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Abstract 

Asymptotic time evolution of a wave packet describing a non-relativistic 
particle incident on a potential barrier is considered, using the Wigner phase- 
space distribution. The distortion of the trasmitted wave packet is determined 
by two time- like parameters, given by the energy derivative of the complex 
transmission amplitude. The result is consistent with various definitions of 
the tunneling time (e.g. the Biittiker-Landauer time, the complex time and 
Wigner 's phase time). The speed-up effect and the negative dispersion are 
discussed, and new experimental implications are considered. 
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Quantum tunneling, where a particle has a chance to pass through classically 
forbidden regions, was one of the first important applications of wave mechanics. 
The total barrier penetration probability may be calculated directly from the sta- 
tionary Schrodinger equation, yet the process time dependence is a more delicate 
phenomenon. The source of the problem is the uncertainty relation: if the incident 
momentum was known exactly, the coordinate would be absolutely uncertain, and 
one could hardly ask a question about the particle transport. In principle, the so- 
lution is evident ffl : a wave packet must be prepared beyond the potential domain, 
far from it, so that one can afford a wide spread, compatible with a relatively small 
momentum uncertainty. After that, one has to wait long enough until the wave 
packet would penetrate through the barrier and leave it completely, being splitted 
in two, forward and backward. However, the problem needs an intricate analysis 
based upon the time-dependent formalism. 

The question: "How much time does the tunneling process take?", has been 
standing since 1932 (e.g. Ref.0). Various definitions, approaches, experiments 
and reviews were worked out (a partial list is given in Refs.|]3||-||22||) in order to 
answer that question which is substantial for physical applications and for a proper 
understanding of quantum theory. 

There are three cardinal approaches to the concept of the tunneling time. 



Path integrals and the semiclassical approximation suggest [H], |20| employ- 
ing a complex time r c . 



A physical clock may be itroduced|23[], at least at the level of a "Gedankenex- 
periment". One way to measure the time spent under the barrier is to study 
the particle spin precession in an external magnetic field j|, which leads to a 
Larmor time r L . Another way is to consider a vibrating barrier || |18[, getting 
the Biittiker - Landauer time r BL . 

One can try to trace the behaviour of the wave packet in interaction with the 
barrier |], UlL Hl> P^l - This approach has two frail points || [7|, ^TJ: i) a 



dependence on the initial state preparation, ii) spreading of the wave packet 
outside the potential region. 

It was found subsequently [13], ^TJ that the definitions of r BL , r L , and r c are 
not inconsistent. On the other hand, for central potential scattering, which is de- 
scribed in terms of partial-wave phase shifts, Wigner^j] (see also in Ref. []!]], Ch. 8) 
introduced a phase time r w (Eq. (1) below) and noted its relation to the causality 
principle. Recent experiments with optical analogs of quantum tunneling|n]-[|T^] 
indicate an importance of the phase time for barrier penetration. 

The purpose of this work is to investigate general features of deformation of wave 
packets in the process of tunneling through potential barriers. It is assumed that the 
initial momentum uncertainty is small, and we look at the large-time asymptotics. 
The result is expressed in terms of the (complex) transmission amplitude A(k), 
which is obtained by the solution of the stationary Schrodinger equation with the 
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energy e = K 2 /2m, the particle mass being m. The total transmission probability is 
given by |A(k)| 2 , and the change of shape of the coordinate probability distribution 
is determined by two time-like parameters related to the energy derivatives of the 
transmission amplitude, 

r w = rf(arg A)/de, r A = d(ln \A\)/de, (1) 

where de = vdn, and v = n/m. (We set % = 1 throughout the paper.) It is 
noteworthy that these parameters appear in other formulae for the tunneling time, 
namely, r c = r w — ir A and r BL = |r c |. The Wigner phase time r w has the 
same form as in the central scattering, where \A\ = 1, in contrast to the problem 
concerned. 

We shall use the Wigner phase-space quasi-distribution, i.e. the Weyl symbol of 
the density matrix (see e.g. a review in Ref.|25|). There are two reasons for that: 
i) this formalism enables one to treat both pure and mixed initial states, and quite 
general types of the wave packets; ii) it is easier to get rid of irrelevant oscillations 
of the wave functions. 

The time evolution of the initial Wigner function po{q,p) can be given by means 



of the phase-space evolution kernel [26] , which represents the fundamental solution 
of the Landau - von Neumann equation for the density matrix. Namely, for any 
time t the Wigner function is 

Pt(q,p) = J L t (q,p;qo,po)po(qo,Po)dqodp . (2) 

We shall consider the scattering problem, where the initial state is prepared with an 
average momentum Pq and has a relatively small momentum dispersion Apo <C \Pq\, 
which are defined as usual, 



P = J ppo(q,p)dqdp, (3) 

(Ap ) 2 = J{p- P ) 2 p (q,p)dqdp. 

Similar definitions hold for the central coordinate Qq and the dispersion Aq . As 
follows from the uncertainty relation, Aq Ap > |, and the inequality may be 
saturated for a set of pure (coherent) states. It is assumed that the potential barrier 
V(q) > is located near the origin and has a finite range D, vanishing for \q\ > D. 
The initial state must be prepared in the free space, which means that \Qq\ — Aq > D. 
The results of the scattering are observed after the wave packets gets out of the 
potential region, say, when t > 2\Q \m/ P . 

In the large time asymptotics, the phase-space evolution kernel for the barrier 
penetration was found [071 to be a sum of two parts, describing transmission and 



reflection, 

L t (q,p;q ,po) x 5(p - p )T(r + ,p ) + 5{p + p )R(r_,p ), (4) 

where r± = qo+tp /m^fq. (If the barrier would be absent, one has T = S(r + ), R = 0, 
and L t is just the solution of the classical Poisson equation for free motion.) The 
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functions T and R have been expressed in terms of integrals involving the transition 
and reflection coefficients, A(k) and B(k), respectively. In particular, the integral 
representation for the transmission propagator is 



T(r,p) = f^Jae-^A{p+ l -a)A{p-^a). (5) 

This representation is manifestly causal. Note that A(k) has the following general 
properties ||28|| : it is analytical in the upper half of the complex plane (having poles 
at Im k < 0), A(k) = A(—R), and lim K ^oo A(k) = 1. Thus the integral in (5) can be 
considered as a contour integral in the complex cr-plane. If q > qo + vt, i.e. q would 
be in advance of the free motion coordinate, the integration contour can be moved 
up to infinity, so T(r,p) = for r < 0. Thus no point of the Wigner distribution is 
transported faster than it would be in the absence of the potential barrier. 

The information transport is somewhat smeared by the fact that the phase space 
distribution is never too local because of the uncertainty relation. Let us consider 
the observable consequences of the causality arguments presented above. We shall 
suppose that the final coordinates are measured, and the detector does not discrimi- 
nate between different momenta. The probability of finding the transmitted particle 
at q, for asymptotically large t, is given by the following integral, 

/oo roo roo ] ] 

dp dr dae- iar A(p + -a)A{-p+-a)p {q-vt + r,p), (6) 
-oo J~oo J—oo Z Z 

where v = p/m. Dealing with this representation, one can make use of the specific 
features of po> which is sharp in p and broad in q. Therefore it is reasonable to 
expand A(±p + \a) in powers of (p — Pq). On the other hand, it was found[^7[ that 
T(r,p) is an exponentially decreasing function of r, so the fact that the r-dependence 
of po is slow may be taken into account. A straightforward calculation shows that 
one can expand the integral in powers of Ap /P . To the first order, the result is 



Here 



V t (q) « \A(P )\ 2 [V°(q) + v r w dV?/dq + v r A 2M t (q)\ . (7) 

V°{q) = J dpp {q - vt,p), M t {q) = J dp{p - P )p {q - vt,p), (8) 

V® represents free motion of the initial wave packet including the usual spreading 
and M t is the first moment of the p-distribution, which also appears in the free 
motion. The time parameters r A and r w in (7) should be calculated by Eq. (1) at 
e = P 2 /2m, and Vq = Po/m. Note that the total transmission probability is 

JdqV t {q)/ JdqV?(q) = \A(P )\ 2 , (9) 

as usual, since the two other terms in (7) are eliminated by the integration. The 
expansion in powers of (Apo/-Fb) can be performed to all orders leading to a sum 
over higher derivatives and higher moments of po with coefficients which depend on 
the barrier shape. 
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Two terms in Eq. (7) describe a distortion of the coordinate propability distribu- 
tion. It is natural to expect that r A > 0, the tunneling probability is increasing with 
energy, so the corresponding term is responsible for an advance in the distribution 
maximum. It is the so-called speed-up effect |TT|| . One can tell that tunneling filters 
out low-energy components of the wave packet. As to the other term, the sign of 
the phase time r w is not definite. (Roughly speaking, it is positive for narrow barri- 
ers and negative for wide barriers). The corresponding term represents interference 
effects, i.e dispersion due to the barrier. This may result in an additional advance 
(if r w < 0), or in a delay (if r w < 0). 

Sometimes the effect of r w is dominating. That is the case when | | is 
constant, as for the example of central scattering, considered by Wignerp4| 



In 



recent experiments with photon wave packets fl4|1, the transmission probability was 



almost insensitive to the wavelength within the light frequency band, so that r <C 
r w . If these parameters are of the same order of magnitude their relative influence 
may be determined by the distance of the detector from the barrier. As the distance 
(i.e. the time t ) is increasing the speed-up effect prevails over the phase-time effect, 
since the slope of the free probability distribution is getting down because of the 
wave-packet spreading. 

Let us consider a simple example where the initial phase-space distribution is 
Gaussian, 

(P-Po) 2 (q-Qo) 



p (q,p) = Cexp 



2(Ap 



2(Ag ) 



(10) 



Here C = (27rA£> Ago) 1 is the normalization constant, and Ap Aq 
is pure. The calculations are straightforward now, 



v°M - 

M t {q) = t 

Q = Qo + 



2nAqC exp 

(g-Q)(A Po ) 

m(Aq) 2 
tv , (Aq) 2 i 



(q - Qf 



2(Aqf 
= (A<?o) 2 - 



| if the state 



(tAp /m) 



In the essential domain, the coordinate distribution given by Eq. (7) is 

V t (q) = \A(P )\ 2 V?(q) [l + v r (q - Q)/{Aqf] , (12) 

where ro = 2tr A (Apo) 2 / m — r w , which may change its sign with time. The max- 
imum of the transmitted distribution is shifted in advance of the that for the free 
propagation by 

AQ = 2v r / ^ 1 + C 2 + l) (13) 

where ( = 2v r /Aq. In the limit of Ap — > 0, Aq — > oo, one has ( 1, and 
AQ w VqTq. Besides, for (Ap ) 2 /m <^ r w /r A t, i.e. if the time is not too large (the 
detector is not too far from the barrier) we get ro ~ — r w and AQ m —vot w in 
agreement with Wigner's prediction and recent experiments jl4}1 . For large positive 
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ro (which may happen at large t), the shift of the maximum is bounded by the 
width of the final distribution. Besides, the transmitted distribution is contracted 
by interaction with the barrier, an effect which is called negative dispersion. 

Equality (7) may be applied to any other distribution, not necessarily Gaussian. 
An example considered in recent microwave simulations of tunneling [11 ] is a step 
function. That distribution cannot be realized in quantum mechanics, but it may be 
considered as a test function. The shift of the half-height point with respect to the 
"free propagation" has been also calculated from Eq. (7). The result is v r h , where 
Th = tr A (Apo) 2 / m ~ rW - Here the speed-up effect is half of that in the Gaussian case, 
Eq. (12). This is consistent with the causality arguments; the wave-packet front 
does not move faster because of a potential barrier. The motion of the peak, as well 
as the motion of the half-height point, have been measured experimentally []14| - fl6| . 

Qualitatively, nonrelativistic particle movement is similar to propagation of light 



signal through a reflecting stack, like in the experiment of the Berkeley group fl4 



The interpretation of the results based upon the light interference picture leads to 



qualitatively similar results p9| . In order to make a quantitative prediction, one has 
to calculate the complex transmission coefficient A(k). If the optical barrier may be 
prepared with a strong frequency dependence, the effect of the "amplitude time" r A 
would be observable, besides the "phase time" r w . Remarkably, the result of the 
measurements would depend on the distance between the detector and the barrier 
region. 

It should be emphasized, in conclusion, that the wave packet shape must be 
taken into the consideration of tunneling through potential barriers. Causality is 
not violated of course, but it manifests itself indirectly in terms of deformation 
of the wave packet. The phase-space formalism, introduced by Wigner, is quite 
appropriate for the investigation, and the main corrections for a wave packet with a 
fairly definite momentum are given in Eq. (7). The resulting effect is an interplay of 
those involving a couple of time-like parameters, defined by Eq. (1) and owing to the 
momentum dependence of the complex transmission amplitude. Thus, the complex 
time r c introduced previously ||, appears actually in the final distribution; its 
real part is coupled to the coordinate derivative of the freely propagating distribution 
and its imaginary part is coupled to the first moment of the momentum distribution. 
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